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Abstract 

By using both thermal field theory and a somewhat more intuitive method, we 
define the electric charge as well as the charge radius of neutrinos propagating inside 
a plasma. We show that electron neutrinos acquire a charge radius of order ~ 6.5 x 
10~^^ cm, regardless of the properties of the medium. Then, we compute the rate 
of plasmon decay which such an electric charge or a charge radius implies. Taking 
into account the relativistic efi^ects of the degenerate electron gas, we compare our 
results to various approximations as well as to recent calculations and determine 
the regimes where the electric charge or the charge radius does mediate the decay 
of plasmons. Finally, we discuss the stellar limits on any anomalous charge radius 
of neutrinos. 
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1 Introduction 



The plasmon decay process is an important cooling mechanism for stars whose core is 
degenerate. It is particularly relevant as regards the evolution of young white dwarves 
and red giant stars. In these two stellar objects, the core is composed of degenerate 
matter, with density p ~ 10^ g cm~^ and temperature in the range from T = 10^ to 10^ 
K. As a consequence of this degeneracy, all electronic levels below the Fermi sphere are 
filled up. A typical value of the electron momentum is pi? ~ 400 keV, so that electrons 
are neither non-relativistic, nor ultra-relativistic. 

The plasmon decay process was first considered a long time ago, in the pioneering work 
by Adams, Ruderman and Woo [|I|, which led to a rich literature 0. The subject has 
recently been revived when it was realized that the plasmon dispersion relations, which 
were used previously, did not incorporate the relativistic behaviour of the electron plasma 
0]. The latest works that deal with these effects are due to Itoh et al. Braaten and 
Segel § and Haft, Raffelt and Weiss §. 

In this paper, we complete the previous analysis and present the calculation of the 
plasmon decay rate in the framework of two different but complementary approaches. 
The first discussion is based on considerations of thermal field theory while the second 
analysis is a more intuitive, kinetic theory-like, method. We show that both reasonings 
lead to the same answer. We also stress the important distinction between the effective 
electric charge which neutrinos exhibit when they couple to a classical electromagnetic 
field and their charge radius which may be strictly defined in the soft regime where the 
energy and momentum of plasmons are small. 

The content of this work is as follows. In section 2, we derive the electric charge 
which a neutrino acquires when it propagates inside matter and interacts with a clas- 
sical electromagnetic field. In section 3, we discuss the charge radius associated to the 
coupling between neutrinos and low-energy plasmons. We comment on these two phe- 
nomenological approximations which lead, a priori, to fairly different descriptions of the 
electromagnetic properties of neutrinos. In section 4, we present the analytic expressions 
for the plasmon decay process. We discuss our numerical results for the case of red-giants 
and white dwarves in section 5 where various approximations to the energy loss rate are 
presented along with their respective domain of validity. We finally conclude that, for 
most astrophysical purposes, the charge radius and the electric charge descriptions of the 
neutrino behaviour, give indeed the same answer. 
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2 The neutrino electric charge and the hard regime 



When neutrinos propagate inside matter, they scatter coherently on charged particles, 
such as the electrons and positrons of the thermal bath. This process is responsible for 
the "additional" inertia which neutrinos possess in a plasma [0. Another possibility is 
that neutrinos diffuse coherently on electrons and positrons which, in turn, couple to the 
electromagnetic field p|. The Feynman graphs responsible for the latter reaction are 
presented in fig. 1 where both electrons (fig. la and lb) and positrons (fig. Ic and Id) 
have been taken into account. Note that the final electronic (or positronic) state is just 
the same as the incoming one, hence a modification of the electromagnetic coupling of 
neutrinos which appears directly at the level of the amplitude of this coherent diffusion. 
The effective coupling of neutrinos to electrons results from the low-energy limit of the 
Standard Glashow-Weinberg-Salam theory : 



C = ^ {^.7a(l-75)^.} {^er{gV-9A75)^e} 



V2 



(2) 



Both Z and W exchanges are taken into account for electron neutrinos so that 

gy = ^ + 2 sin^ 9w and = ^ , 

whilst, for z/^ and Ur where the Z boson alone mediates the interaction, the coupling 
parameters are 

2 sin^ 9w and gA = — - . (3) 



9v 
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-- + 2 sm Ow and g^ = -- . 

In order to compute the effective coupling of neutrinos to the electromagnetic field 
inside a plasma, the spectator electrons must be integrated out. Let us consider a vol- 
ume V of matter. The electronic state \E'e describes the presence of a single electron of 
momentum p in the whole volume V, hence the normalisation 



m 



VE 



u{p, 



(4) 



where m and E stand respectively for the mass and the energy of the electron. For 
positrons, u{p , s) is replaced by v{p, s). The electromagnetic coupling of neutrinos may 
therefore be expressed as a sum of the amplitude of the coherent diffusion pictured in 
fig. 1, over the electronic and positronic states of the thermal bath : 



Ve{p) 



^ {^.7.(1-75)^.} A, J 
^e(p , s)7" {gv - 9Al5' 
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Veip) 



^e-(p , S)e7^— — : 



-l" {gv - 9Al5) ^e(p , S) 
1 



Here, = {E,p) denotes the electron (or positron) momentum whilst = [uj, q) refers 
to the plasmon. Since the initial and final states of the spectator electrons are the same, 
expression (^) may be written as a sum of traces. The statistical occupation numbers of 
electrons and positrons are respectively taken care of by the Fermi-Dirac functions rie{p) 
and r]e{p). The effective Lagrangian (|) becomes 

eG 



C 



eff 



Veip) 



{^,7, (1 - 75) Ap J dp X 

ip + m) 7- [gv - gAl^) [p + Q + m] 7^ 



(6) 
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Q2 + 2P-Q 



(Z' + m) 7^ (/^ - ^ + m) 7" {gv - gAl^) 



Tr 



g2 - 2P-Q 

iP - m) 7^ + m) 7" {gv - gAl,] 



- Ve{p) ' 



+ 



Tr 



Q2 + 2P-Q 

iP - m) 7" {gv - 9A15) {-P + Q + m) 7^3 



- 2P-Q 

where dp is the Lorentz-invariant differential element 

After the reduction of the above-mentionned traces and some straightforward algebra, the 
effective coupling of neutrinos to the electromagnetic field may be expressed as 

^e// = ^ {^.7.(1-75)*.} r"^ A, , (8) 

where the tensor F"^ may be separated into its symmetric and antisymmetric components : 
T^^ = Aegv [ dp Mp) + ri,{p)} 



{P-Q) (P'^Q^ + P^Q'') - {P-Qf g''^ - Q'^P'^P^ 



2tegAe''''^'' Jdp{Ve{p) - Ve{p)} 



{P-Qf - {Qy2f 



{P-Qf - {Qy2f 



(9) 



This effective coupling F"^ has been derived fairly intuitively. As pictured in fig. 1, the 
neutrino wave interacts coherently with the charged species of the plasma - for that matter 
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mostly electrons and positrons since the contribution of nuclei is negligible. Expression 
(^) may be directly obtained from mere considerations of thermal field theory. 

At the one-loop level, the Feynman rules for finite temperature field theory are quite 
simple 0. The only ingredient is here the electron propagator, given by 



Sf{P) = ip + m) 



.P2 _ 772^ + ie 

where the statistical factor Nf{P) may be expressed as 



2n6{P^ - m^)NF{P) 



(10) 



N (P) = , ^(--Po) .11) 

eiPo - /^)/T + 1 + ,{-Po + /x)/T + 1 • 

Using the same conventions as before, the only non- vanishing contribution to the one-loop 
diagram of fig. 2 is 

= e| ^ Tr [7" ((7y - ^7^75) (yP + m) 7^/^ + ^ + m) 
j 6{P'-m')NF{P) 6[{P + Qf-m']NF{P + Q) \ 

\ (P + g)2-m2 P2_^2 J ' ^ ) 

which, after a change of variables, reduces to relation (^. The same result could also 
have been derived in the imaginary-time formalism. 

The effective electric charge of the neutrinos which propagate inside a plasma may 
be defined as the limit of expression (^) where describes a classical electromagnetic 
field. The behaviour of finite temperature field theory and especially of its perturbative 
approach of interactions is quite rich. Inside matter, photons no longer propagate at the 
speed of light. On the contrary, they follow dispersion relations which imply the exis- 



tence of longitudinal and transverse modes - see equations (^) and (0). The transverse 
photons are extremely similar to the ordinary quanta of the vacuum. They should be 
disentangled from the longitudinal modes which are mere collective excitations of the 
plasma. That is why the transverse branch alone exhibits an asymptote where the energy 
uj and the momentum q of plasmons reach up to infinity at constant = lo^ — . In 
this so-called hard limit |[T^, the transverse photons behave as if they propagated in the 
vacuum, with the non- vanishing mass my. Note that is always of order cJq, where 
the plasma frequency cuq is smaller than the electron mean energy < E'e > by a factor 
of e = A/ivra, i.e., cuq ~ e < -Eg >~ eT or fpp. Therefore, in the hard regime, the 
components uj and q may vary regardless of which behaves as a small perturbation of 
order a. In relation (||), expressions such as (P-Q) or P^Q^^ may be considered as leading 
contributions with respect to the correction ~ 0(a). The tensor F"^ may therefore 
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be expanded as a function of the small parameter Q^, and the leading term turns out to 
be the symmetric expression 

T " Aegv J dp{rje[p) + Ve{p)} | | 

+ terms of order a . (13) 



As expression ([T3| ) is folded into the vertex (H), it further simplifies. In the Lorentz gauge 
where dpA^ = 0, the contribution Q^Ajs vanishes. Note that the effective interaction 
(I) is gauge invariant since F^^Qq, = F^^Q/? = 0, as is obvious from expression (H). 
Therefore, restricting ourselves to the Lorentz gauge does not alter the general import 
of our reasoning. Finally, neutrinos may be considered as massless species, a fairly good 
approximation since plasmon decay merely involves the electron neutrino for which gv — 1 
and whose mass is orders of magnitude smaller than the typical energies at stellar cores. 
Since expression ^i^^ (1 — 75) "^i, vanishes in that limit, the neutrino coupling further 
simplifies into the effective Lagrangian : 

Ceff = - eg. {^.^^L^.} Ap . (14) 

The left-helicity operator L is defined as L = (1 — 75) /2. The neutrino behaves as if it 
had the electric charge 

Q, = AV2GFgv^p . (15) 
This charge depends on the properties of the plasma through the integral 

In a degenerate gas of electrons, a typical situation inside red-giants or white dwarves, 
this plasma integral over the statistical distribution of electrons may be expressed as 



where vp = Pf/Ef denotes the velocity at the surface of the Fermi sea. At the helium 
core of red-giant stars, for instance, the density reaches up to ~ 10^ gcm^'^. The Fermi 
momentum is pf ~ 410 keV, hence an effective electric charge ~ 6.4 x 10^^^ which, 
in spite of being some 13 orders of magnitude smaller than the charge of the electron, 
induces a significant energy loss, so large that it dominates the evolution of stars during 
their ascension of the red-giant branch. 
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3 The neutrino charge radius and the soft hmit 



In the low energy part of the transverse branch as well as all along the longitudinal branch, 
the energy uj and the momentum q of plasmons are much smaller than the average energy 
of the charged species of the thermal bath. In this so-called soft regime where q is 



small, the components of the plasmon four-momentum may be considered of order 
e = v47ra with respect to the typical electron momentum, i.e., 

1^ ~ C(e) ~ 0{V^) < 1 . (18) 

Notice that in the hard regime, is of order 0{a) whereas both components and 
P'^ are leading order terms. On the contrary, in the soft regime which we discuss now, 
expressions such as q or y/oj'^ — q^ are of order C(e), so that the perturbative expansion 
of expression (|^) is now quite different than in the hard-loop approximation. In particular, 
the axial part of the tensor F"^ is of order 0{e). However, it does not interfere with the 
vector component of F"'^ as regards plasmon decay, and therefore contributes a term only 
of order 0{q) which we will not consider in what follows. The leading contribution to 
F"^ is the vector component : 



F"^ = ^egv JdpMp) + r],ip)} 



(P-Q) (P'^Q^ + P^Q°) - {P-Qf g""^ - Q^P^'P'^ 



{P-QY 

+ terms of order e . (19) 

As noticed already some thirty years ago |jl| , the previous expression is merely the polar- 
ization tensor 11°"^ associated to the propagation of electromagnetic waves inside a plasma, 
i.e., 

_ ^n"'^ , (20) 

e 

hence the effective Lagrangian 

C.jf = ^^^^ {^.7.L^.} A, , (21) 

where the electromagnetic field couples to the neutrino current through the polarization 
tensor. Note that 11°^ may be broken up into the operators P^^ and P^^ which respec- 
tively project the polarization on the transverse and longitudinal modes |Tl| : 



U^P = UT{u;,q)Pf + UL{uJ,q)Pf . (22) 
Since the plasmon four-momentum follows the dispersion relation 

Q2 + nT,L(a;,g) = , (23) 
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expression (pTf) simplifies further into : 



C.ff = - ^^^^ {^ul'L^.} A, . (24) 

Tlierefore, in tlie soft regime, neutrinos interact with the electromagnetic field as if they 
had a non-vanishing charge radius : 

2 

Ceff = -e^-^Q^ [^>,^fm>J] , (25) 

with the plasmon four-momentum factored out in the coupling, as is typical. This 
thermal charge radius, which neutrinos built from their coherent diffusion on the spectator 
electrons of the plasma, may be expressed as 

< > = 9v , (26) 

and may be understood as the superimposition of two distributions of opposite electric 
charge, with same centers, but with different spatial extensions. By surfing on the elec- 
trons of the plasma, the neutrino behaves as if it acquired an inner structure whose typical 
size is surprisingly independent of the properties of the thermal bath : 



(6.5 X 10"^'' cm)^ gy . (27) 



We therefore conclude that, in the soft limit, neutrinos interact with the electromag- 
netic field through a charge radius coupling whereas, in the hard regime, they possess an 
effective electric charge. 

In a degenerate medium, the plasma frequency ioq is related to the Fermi momentum 

Vf 



= ^ ^ . (281 



4a 'Pp 



■^0 



Stt Ep 

A fully relativistic treatment of the electron bath leads to the dispersion relation for the 
longitudinal mode |l2l : 

where vp = pp/Ep stands for the velocity at the surface of the Fermi sea. The energy 
and the momentum of plasmons are respectively denoted hj ujl and q. The longitudinal 
branch extends up to the point where it crosses the light cone, for 



4f- >" 

Vp \2vp 



l+Vp 
1 — Vp 



1/2 

- 1 ) ^ uJo . (30) 
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As q increases from to qmax-, the four-momentum steadily decreases from the plasma 
frequency ujq^ down to 0. Note that and q cannot reach up to infinity, except in the 
ultra-relativistic case where, nevertheless, they cannot vary regardless of the value of Q^. 
The longitudinal branch corresponds therefore to the soft regime. In particular, when 
electrons are non-relativistic, i.e., for vanishing vp, relation (|29|) simplifies into 

u}l = Uq , (31) 

and the longitudinal branch extends up to qmax = ^o- The non-relativistic regime is a 
good example of soft behaviour in so far as ujl and g, which cannot exceed ujq, are obviously 
of order 0{e) with respect to the Fermi momentum pp- As regards the transverse mode, 
the dispersion relation inside degenerate matter may be expressed directly as a function 
of the effective plasmon mass : 



The low-energy part of the transverse branch corresponds to the soft regime, whereas the 
hard limit is recovered for cuj- ^ g — +oo. When the hard regime obtains, transverse 
plasmons behave as if they propagated with the effective mass 

As q increases along the transverse branch, note that the four-momentum raises from 
uoq^ to m\. As is clear from relations (0), (|28|) and (0), this asymptotic mass may be 
readily expressed as a function of the integral 1p which dominates the behaviour of the 
neutrino electric charge Qy : 

m| = Ae'^Xp . (34) 

When the stellar plasma is ultra-relativistic, the transverse mass of high-energy {i.e., 
hard) plasmons is rriT = ^^3/2 Uq. On the contrary, for non-relativistic electrons, relation 
(P^) becomes 

cut^ = cJo^ + , (35) 

and plasmons propagate with the mass ujq all along the transverse branch, even in its soft 
region. 

Both longitudinal and transverse dispersion relations were first derived by Jancovici 
13| , but had a much more complicated form. It was realized later that the simplified 



expressions (^) and (|32D could be used in a much wider regime [§, |T2 



A thorough inspection of relation (^ shows that equation (|T3p behaves as the mere 
limit of equation (|T9|) for vanishing Q^, at fixed plasmon momentum. The transition 
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between the soft and hard behaviours is therefore smooth. The expressions derived in 
this section are surprisingly more general than what could be naively guessed. Not only 
are they valid in the soft limit, but they cover the entire plasmon spectrum and, as a 
consequence, the hard regime as well. The effective charge radius is a good description 
of the electromagnetic properties of neutrinos. To illustrate this point, let us consider 
the effective interaction (plf) which has been strictly derived in the soft regime. For 
high energy plasmons of the transverse mode, may be replaced by 4 Xp so that the 
charge radius coupling (< > Q^/6) translates into the electric charge of relation 
(|T5|). The charge radius, which has been defined in the soft regime, pleasantly behaves 
as an effective electric charge when the hard limit is taken. More generally, it translates, 
along the asymptote of the transverse branch, into the electric charge 

Gpgv 2 /Q«\ 
= — ■ (36) 

The Lagrangian (|2^) provides therefore a comprehensive description of the electromag- 
netic properties of neutrinos propagating inside a plasma, valid for the entire spectrum 
of both longitudinal and transverse modes. It reduces to an electric charge when the 
hard regime applies, ie., in the asymptotic limit of the transverse branch where plasmons 
behave classical field. 



4 The plasmon decay rate 

Inside stars, plasma waves may decay into neutrino pairs as a result of the effective 
coupling (^). At the helium cores of red-giants where the stellar density is fairly large, 
the plasma frequency reaches up to cuq ~ 20 keV. Once produced, the neutrinos do not 
interact with stellar interiors. Most of the stars are indeed completely transparent to 
these species which, therefore, escape freely from their cores. The corresponding cooling 
rate depends on the distribution of plasma waves as well as on the efficiency with which 
a plasmon is transmuted into a neutrino pair. 

The decay rate of plasmons obtains from the phase space integral 

= ^ / dK d\ (27r)^ 6' {Q~K~ h) Y: \M (7Q ^ lyiy) T , (37) 

where the relevant matrix element is 

Mi^Q^uu) = ^/2^£^g2 ^u^k^yiVih)} . (38) 

In the Lorentz gauge, the polarization of transverse plasmons are the pure vectors er, 
transverse to the direction q of motion, whilst the longitudinal polarization is taken care 
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of by the four- vector 



= {q.-^^L)/^|u:l-q^ . (39) 



A straightforward calculation leads to the decay rate 

1 GWv 1 



D 



h m ■ («) 



487r2 a u 

In the hard regime, where the electromagnetic properties of neutrinos are well described 
by the electric charge eQy - see relations ([l^) and (^) - the decay rate becomes 



The neutrino energy loss may be derived from the sum over the Bose-Einstein distribution 
risiuj) of the individual decays 

(Pq 



= J^9h Vb{uj) uZVd . (42) 

Each neutrino pair carries outside the star the energy uj of its progenitor. The helicity 
factor Qh respectively takes the values 1 and 2 for the longitudinal and transverse modes. 
The additional factor Z corrects for the term l/2cij in relation (P?!). For ordinary particles 
with mass m, the relation between the energy E and the momentum p is merely 

- p 2 ^ p2 ^ , (43) 

so that, when the Lorentz invariant expression d^P 5 {P^ — m^) is integrated over the 
positive values of P°, it translates into the conventional phase space element d^p / 2E. 
For plasmons, expression (^3]) is replaced by the more intricate dispersion relations (|23|) 
which, once integrated over the positive values of become 

d'Q e{Q') 6{q' + n(Q°,g)} = Z |^ , (44) 



with 

z-' = ' " 



For the transverse mode of a degenerate plasma, this normalisation factor may be ex- 
pressed as 

n 2 ( 2 2 2 



3 ojq^ojt^ + {ojt^ + g^) {ujt^ — f ^g^^ — 2ujt^ {ujt^ — q' 
whereas, for the longitudinal branch, it becomes 



(46) 



2cjr2 I [uJl^ - vlq^ 



Q' 13^0^ - {u,' -vW) 



(47) 
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In the non-relativistic regime for which the velocity vp ^ 0, Zt merely reduces to unity 
whereas Zl tends towards ujq^ /Q"^. 

Once again, similar expressions for the energy loss may be obtained in the thermal 
field theory approach. The neutrino (antineutrino) production rate, due for instance to 
transverse photon decay, is easily derived by using the cutting rules of Kobes and Semenoff 



14| , given by 
dN, 



X 27r5 

Here K and Q are respectively the neutrino and photon four-momenta. The transverse 
photon projection operator is 

Vf^ = - + q'q^/q^ , (49) 



with all its other components set equal to zero ||rT|. After performing the integrations, 
the result takes the final form 

R, = ^rt^Zr,ui^)(Qi . (50) 
The energy loss rate may be inferred directly from the previous integral : 

= ^ r^'dqZrVBico) (gf , (51) 
487r^a Jo ^ ^ 

where the factor F^, takes into account the three neutrino species produced in plasmon 

decay 

= Y.^gU^) = ^-2 sin2^vK + 12 sin^^iy • (52) 
As regards the longitudinal branch, a similar expression obtains where the plasmon mo- 
mentum is integrated up to qmax 

In the next section, the behaviour of the exact cooling rate e = + will be discussed 
in various stellar environments. Two approximations of this exact energy loss may be 
defined : 

(i) The estimate obtains when the charge radius (^) is replaced by the electric 
charge (p!^), so that the transverse cooling is now 

Q IQaG^F^ ^2 



2 



poo 

ll / q^dq Zt VBi^j) , (54) 
Jo 
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whereas the longitudinal energy loss is 



4-'-^il.l^""i-<liZ.mMQ' . (55) 

(ii) The non-relativistic approximation corresponds to the limit where vp vanishes. 
The plasma frequency is therefore defined as the ratio 

2 4a pI 

UJn = . 56 

As the dispersion relations (^TJ) and (P^D now apply, the expressions for the energy loss 
rate considerably simplify. We simply recover the results of Adams, Ruderman and 
Woo later on corrected by Zaidi 0. Neutrino emission by transverse photons occurs 
at a pace : 

where the photon density is defined by 

1 ^ 



vr^ J 



q dq tib^ujt) ■ (5^ 



As regards the energy dragged away from the longitudinal mode, the corresponding loss 
rate may be expressed exactly 

= tII^ Vsiu^o) ■ (59) 
IzoU TT^a 



5 Discussion 

Neutrino cooling dominates the slow contraction of red-giant stars |jl5|, whose core tem- 



perature steadily increases in time, up to the point where, at the tip of the red-giant branch 
of the Hertzsprung-Russell diagram, the triple-alpha reaction sets on. For low-mass stars, 
typical conditions at the center are p ~ 10^ g cm~^ with T ~ 10^ K. During this stage of 
their evolution, stars have a burning shell of hydrogen surrounding a degenerate helium 
core whose mass increases at a pace dominated by the neutrino losses. Thus, the mass of 
helium processed on the red-giant branch may be related to plasmon decay and its asso- 
ciated cooling []TB|. As a matter of fact, observations of the color- magnitude diagrams of 



globular clusters offer an interesting handle on the cooling of red-giant cores by neutrino 
emission. The latter mechanism also dominates the early evolution of white-dwarves, 
slightly after their formation, when they are still relatively hot, with central temperature 
T ~ 10^ to 10^ K. The cooling of white-dwarves is expected to be measured in the near 
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future, using asteroseismology techniques. The slow variation of the period of pulsating 



white-dwarves may be directly related to the decrease of their central temperature [17 



The case of red-giant stars is embraced in fig. 3a where the exact energy loss rate e 
varies with the density p, at fixed temperature T = 10^ K and fixed electron fraction 
Ye = 0.5 (pure helium). For small values of the density p, the Fermi momentum 

PF = {3rf^Yep/muy^' (60) 

is smaller than the electron mass m, and the non-relativistic regime applies. The point 
p = 10^ g cm~^, for instance, corresponds topp = 190 keV, so that the plasma frequency is 
ti^o — 6 keV. Even at such a low density, electrons are degenerate because the temperature 
T ^ 9 keV is negligible with respect to the Fermi momentum pp. Since T exceeds the 
plasma frequency uq in the left-hand side region of fig. 3a, relations (|57|) and (|^) simplify 
respectively into 

e^T > coo) = ^ GlF^ , (61) 

and 

e^«(T > uo) = T . (62) 

^ ^ 1260 vr^a " ^ ' 

The transverse rate overshadows completely the longitudinal emission. Note that at low 

density, the energy loss e ^ e^^ + e^^ increases with p. On the contrary, in the high 

density regime (right-hand side domain of fig. 3a), electrons are ultra-relativistic and 

the temperature is now much smaller than the plasma frequency, hence a decrease of e 

with p. The Boltzmann factor exp(— c^o/T) drops as there are fewer and fewer plasmons. 

The maximum which the curve exhibits corresponds to the transition between these two 

regimes. 

In fig. 4a, devoted to white-dwarf cooling, the density has been set equal to the typical 
value of 1.8 tons cm~^, a pure carbon composition has been assumed (Yg = 0.5), and 
the temperature is varied from 10^ up to 10^ K. The Fermi momentum pp ~ 500 keV 
corresponds to the relativistic velocity vp = 0.7 and to the plasma frequency luq = 23.2 
keV. Note also that the transverse mass of high-energy photons is = 24.7 keV. Since 
the maximal temperature considered here, T = 86 keV, is much smaller than the Fermi 
momentum, the gas is still degenerate. The increase of the energy loss rate e with T is 
impressive. 

The conditions of density and temperature in the plots 3b and 4b are respectively the 
same as in fig. 3a and 4a. The solid line stands for the ratio e^/e and corresponds to the 
electric charge approximation. The short dash curve refers to the non-relativistic estimate 
e^^/e. The numerical fit published by Itoh et al. is extensively used for astrophysical 
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purposes. Its ratio to the exact cooling rate is represented here by the dotted curve 
while the short dash-dotted line stands for the recent numerical estimate by Haft et al. 
0]. The energy loss rate e is a function of both the density p and the temperature 
T. In a degenerate medium, which is the situation under scrutiny here, the density alone 
determines the Fermi momentum, the plasma frequency as well as the dispersion relations 
of the plasmon propagation. Therefore, the structure of both transverse and longitudinal 
modes depends only on the parameter p. As for the temperature, it conditions the filling 
of the plasma states along the various branches and governs, for instance, the above- 
mentionned density n^. Some remarks are in order. 

(A) In the regime where the temperature exceeds the plasma frequency, the high- 
energy tail of both transverse and longitudinal branches are filled up. Note that ujl^ — 
vanishes near the tip of the longitudinal branch. Therefore, the decays of longitudinal 
plasmons are suppressed and neutrino cooling is dominated by the behaviour of the high- 
energy, i.e., hard, transverse photons. The loss rate e*^ is a good estimate of e in the 
high temperature regime. That is why the ratio e*5/e (solid curve) becomes unity in the 
left-hand side of fig. 3b and in the right-hand portion of fig. 4b. 

(B) At low temperature, i.e., in the limit where ujq exceeds T, the low-energy plasma 
states only are thermally excited, and the average momentum q is, a priori, much smaller 
than the plasma frequency. Plasmons may be pictured as if they were at rest, with energy 
ujq, vanishing momentum q and therefore an effective mass of order the plasma frequency 
ujq. The exact decay rate ( ^DD varies as (Q^)^ whereas the approximate relation (|T) is 
proportional to m^Q^. Therefore, when T is small with respect to cuq, the ratio e*^/e 
behaves roughly as (mr/cjo)^. In fig. 3b for instance, the sohd curve tends asymptotically 
towards the ultra-relativistic limit of (3/2)^ = 2.25 when p is large. In the low temperature 
part of fig. 4b, it converges towards the value (24.7/23.2)^ ~ 1.3 (see above). There is 
nevertheless an exception to this behaviour. When the plasma is non-relativistic, the 
longitudinal branch is noticeably flat and, unlike for the other mode, all its propagation 
states are equally filled up, even at extremely low temperature. The momentum q is 
not forced to vanish because the Bose-Einstein factor rjBi^^o) is constant throughout the 
longitudinal branch. In this regime, the energy loss 

e^T < -o) = ujo' e-^o/T ^gg) 

1260 TT^a 

eclipses completely the transverse production of neutrinos 

e^T < -o) = 0..'"' T^'' e-^o/T . (54) 

48v2 vr'/^ a 

Therefore, at low temperature and for a non-relativistic degenerate electron gas, the 
electric charge approximation overestimates the exact rate by a factor of 35/8 ~ 4.4. 
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(C) At low density, the non-relativistic regime obtains and is a good estimate of 
the exact loss rate. Consequently, in the left-hand side of fig. 3b, the short dash curve 
approaches unity. On the contrary, for large values of the density, the non-relativistic 
expression (^) overestimates slightly the plasma frequency uq. In this regime, the loss 
rate e behaves as exp(— cuq/T"), hence the exponential drop which the ratio e^^/e exhibits 
as the density increases (right-hand side of fig. 3b) or as the temperature diminishes 
(left-hand side of fig. 4b). 

(D) For red-giant stars, the fit which Haft et al. propose is impressively good. Notice 
how the short dash-dotted curve remarkably sticks to unity throughout fig. 3b. The 
Japanese approximation slightly underestimates the exact rate below p = 10® g cm~^. As 
regards white-dwarf cooling. Haft's fit is still fairly good, especially between 10^ and 10^ K. 
On the contrary, the result by Itoh et al. exhibits wiggles at the transition temperatures 
of the numerical estimate whose quality is not completely satisfactory. Furthermore, 
below T = 10'' K, the corresponding dotted curve diverges whereas Haft's fit still gives 
an accurate prediction. 

The coherent diffusion of neutrinos on the electrons of a plasma induces an effective 
interaction with the electromagnetic field. An intuitive approach as well as a more formal 
discussion performed in the framework of thermal field theory have been presented here. 
Both methods lead to the same conclusions. Neutrinos which propagate inside matter 
acquire an effective charge radius of order 6.5 x 10~^® cm Qy'^, regardless of the specific 
properties of the plasma. The description of the electromagnetic interactions of neutrinos 
by a charge radius, strictly valid in the limit of the soft-loop regime, turns out to be more 
general and may apply to the entire spectrum of plasmons. In particular, in the hard 
limit of high-energy transverse quanta, this neutrino charge radius translates into a mere 
electric charge. Since plasmon decay is an efficient mechanism of stellar energy loss when 
the cores of stars are hot, i.e., in a regime where the estimate e*^ is accurate, the thermal 
electric charge Q^, of neutrinos is therefore relevant to astrophysics. However minute as it 
may look, it is nevertheless responsible for the evolution of red-giant stars as well as for 
the early cooling of nascent white-dwarves. 

Suppose finally that the tau neutrino has a large charge radius of order 6.5 x 10~^® 
cm. Plasmon decay into u^- pairs would induce an extra neutrino emission comparable 
to the standard production. Observations of the magnitude shift between the tip of the 
red-giant branch and the horizontal sequence of the HR diagram preclude the existence 
of any additional neutrino cooling with a magnitude larger than the standard energy loss. 
We readily infer a limit of order 6 — 7 x 10~^® cm on any anomalous charge radius. For 
electron and muon neutrinos, the laboratory constraints, of order 1 — 3 x 10~^® cm, are 
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already more stringent |T^. However, for tau neutrinos, the stellar limit compares fairly 
well to the upper bound of 8 — 9 x 10^^^ cm derived from the ASP and CELLO experiments 
0. 
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Figure Caption 



Fig. 1 The neutrino current diffuses colierently on the electrons (la and lb), and on 
the positrons (Ic and Id) of the thermal bath. Through its interactions with the 
spectator particles, this current couples to the electromagnetic field A^, hence an 
effective electric charge (or charge radius) for the neutrinos which propagate inside 
a plasma. 

Fig. 2 One-loop Feynman diagram contribution to the electromagnetic interactions of 
neutrinos, in finite temperature field theory. 

Fig. 3 The total neutrino emissivity e as a function of the plasma density p for pure 
helium, at T = 10^ K. Fig. 3a shows the exact absolute rate whilst fig. 3b presents 
the ratios between various approximations and the exact rate. 

Fig. 4 The neutrino emissivity as a function of the plasma temperature for pure carbon, 
at p = 1.8 X 10^ g/cm^. Fig. 4a shows the exact absolute rate whilst fig. 4b exhibits 
the ratios between various approximations and the exact rate. 
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